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We generalize the kinetic theory of fluids, in which the description of fluids is based 
on the geodesic motion of particles, to spacetimes modeled by Finsler geometry. Our 
results show that Finsler spacetimes are a suitable background for fluid dynamics and 
that the equation of motion for a collisionless fluid is given by the Liouville equation, 
as it is also the case for a metric background geometry. We finally apply this model 
to collisionless dust and a general fluid with cosmological symmetry and derive the 
corresponding equations of motion. It turns out that the equation of motion for a dust 
fluid is a simple generalization of the well-known Euler equations. 
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1. Motivation 

Modern cosmology has confronted the physics community with a number of un¬ 
explained observations, such as the accelerating expansion of the universe and the 
isotropy of the cosmological microwave background. The standard model of cosmol¬ 
ogy aims to explain these observations by the presence of dark energy in the ACDM 
model and an inflationary phase in the early universe. However, the mechanisms 
behind these explanations, i.e., the constituents of dark energy and the driving force 
of the inflation, are not yet understood. This opens the possibility to create various 
models both from particle physics and gravity theory. The model we discuss here 
belongs to the gravitational category. The basic idea of this model is to replace the 
metric geometry of spacetime from general relativity by a Finsler length measure. 
In this article we derive a theory of fluid dynamics based on this Finsler geometric 
background. 

The work we present in this article has two main ingredients. The first ingr edient, 
which serves as the background geometry, is the Finsler spacetime framework.tHSEIl] 
This framework introduces a notion of Finsler geometry which provides a unified 
description of a Lorentzian causality, observers and gravity, and which can be used 
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as a background for field the ories such as electrodynamics. The second ingredient 
is the kinetic theory of fluids This theory is based on the idea that fluids can 
be modeled as being constituted by particles, whose worldlines are geodesics on 
the background spacetime. In the continuum limit this theory yields equations of 
motion for a density function on a subspace of the tangent bundle of spacetime. In 
this article we generalize this kinetic theory of fluids to Finsler sp acetimes. I n this 
model the density function becomes a function on observer space,^ ^ l ^^ l ^^ l which 
is the space of physically allowed four-velocities. We apply this formalism to two 
physically motivated examples. 

The outline of this article is as follows. In section [2] we provide a brief introduc¬ 
tion to the geometry of Finsler spacetimes. From this starting point we construct the 
space of physical observer velocities in section |3] and discuss its geometric structure. 
We then use this structure to construct a model of fluid dynamics based on kinetic 
theory in section |4l In section [5] we apply this model to two particular examples: 
a collisionless dust fluid and a fluid with cosmological symmetry. We end with a 
conclusion in section m 

2. Finsler spacetime geometry 

We start our discussion with a brief review of the background geometry. In this 
section we introduce the basic geometric objects on a Finsler spacetime, which 
will later be used f or the construction of fluid dynamics. Our starting point is the 
following definitionPSEl 

A Finsler spacetime (M, L, F) is a four-dimensional, connected, Hausdorff, para- 
compact, smooth manifold M equipped with continuous real functions L,F on the 
tangent bundle TM which has the following properties: 

(1) L is smooth on the tangent bundle without the zero section TM \ {0}. 

(2) L is positively homogeneous of real degree n > 2 with respect to the fiber 
coordinates of TM, 


L{x,\y) = L{x,y) VA > 0 


( 1 ) 


and defines the Finsler function F via F{x,y) = \L{x,y)\^. 

(3) L is reversible: \L{x,—y)\ = \L[x,y)\. 

(4) The Hessian 



( 2 ) 


of L with respect to the fiber coordinates is non-degenerate on TM \ X, 
where X C TM has measure zero and does not contain the null set {{x,y) € 
TM\Lix,y) = 0}. 
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(5) The unit timelike condition holds, i.e., for all a; G M the set 

ria; = jy G T^M \L{x,y)\ = 1, 

r , , , , . L(x, y) 1 , . 

9abixi y) has signature (e, -e, -e, -e , e = —--^ ) 3 

\L{x,y)\ J 

contains a non-empty closed connected component Sx Q 0,^ G T^M. 

Here we have used coordinates (a;“) on M and the induced coordinates (x°',y°') 
on TM such that {x,y) = y°‘da\x G TxM. The Finsler function F introduced above 
defines a length functional 

s[y] = J dr F{j{T),j{T)), (4) 

which measures the length of a curve t i—>■ 7 (t) on M. An important class of curves 
called Finsler geodesics are those for which the length functional becomes extremal. 
They satisfy the geodesic equation 

r + N\{j,jN = 0. (5) 

Here we have introduced the Cartan non-linear connection 

= Idb - d,F^)] , ( 6 ) 

where da = d/dx°', da = djdy^" and g^ °-b is the inverse of the Finsler metric 

9ab{x,y) = \dadbF^{x,y). (7) 

The Cartan non-linear connection induces a unique split of the tangent bundle 
TTM into horizontal and vertical parts, TTM = HTM © VTM. The horizontal 
tangent bundle FITM is spanned by the vector fields 

{<1, = a, - N\di ,}, (8) 

while the vertical tangent bundle VTM is spanned by {9o}- The corresponding 
basis {do, 9a} that respects this split is called the Berwald basis. Its dual basis of 
T*TM is given by 

= dj/“ + iV“f,dx^}. (9) 

Using the Berwald basis we can find another description for Finsler geodesics. For 
this purpose we consider the canonical lift F : R —>■ TM of a geodesic 7 to the tan¬ 
gent bundle TM. Writing F in coordinates (F“,r“) = ( 7 “, 7 “) we find the geodesic 
equation 

f“ = f“ , f“ = -iV“b(r“, f“)f''. (10) 

Since this is a first order differential equation, the canonical lifts of Finsler geodesics 
take the form of integral curves of a vector field S = y°‘Sa on TM, which we call 
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the geodesic spray. Finally, the Berwald basis and the Finsler metric allow the 
construction of a metric 

G = -gf, 0 6y^ . (11) 

on the tangent bundle, which is called the Sasaki metric. 


3. Observer space geometry 

After discussing the most important Finsler geometric structures on the tangent 
bundle TM in the previous section we now restrict ourselves to a particular subspace 
of TM. Recall from the definition of a Finsler spacetime that for each x G M there 
exists a shell Sx C TxM of future unit timelike vectors, which corresponds to the 
four-velocities of test masses and physical observers. Their union 

O = U (12) 

xGM 


is called observer space. This is the space on which we will define fluid dynamics in 
the remainder of this work. 

The observer space is a seven-dimensional submanifold of the tangent bundle 
TM and thus inherits a number of geometric structures from the Finsler geometry. 
In particular, the geodesic spray S introduced in the previous section is tangent 
to the observer space, and thus restricts to a vector field r on O, called the Reeh 
vector field. This relation between the geodesic spray and the observer space is in 
fact necessary for the interpretation of O as the space of physical four-velocities: it 
means that a test mass possessing a physical initial four-velocity and following a 
geodesic, and thus an integral curve of S, will have a physical four-velocity at all 
times. The canonical lifts of all physical geodesics, corresponding to freely falling 
test masses, are thus given by integral curves of the Reeb vector field r on O. 

The second important structure on the observer space is the restriction G of 
the Sasaki metric G to O. This also equips O with a volume form E = VoIq. This 
volume form has a number of properties which are relevant for the construction of 
a theory of fluid dynamics. The most relevant property for the remainder of this 
article is the fact that it is preserved by the flow of the Reeb vector field, CfiT = 0. 

From the Reeb vector field r and the volume form E we finally define the particle 
measure w = irS. It is the unique six-form (up to a constant factor) which has 
the following properties necessary for the construction of fluid dynamics from the 
kinetic theory. Most importantly, its restriction to a hypersurface a C O which is 
not tangent to the Reeb vector field r is a non-vanishing volume form on a. Further, 
it is closed, doj = 0, and preserved by the flow of the Reeb vector field, £rW = 0. 
The relevance of these properties will become clear in the following section, when 
we use it as an ingredient to develop a theory of fluid dynamics. 
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4. Kinetic theory and fluid dynamics 


We now turn our attention to the kinetic theory of fluidJ^^^^ and apply it to fluids 


on the Finsler spacetime background geometry detailed in the preceding section. 
For this purpose we first briefly review how fluids can be modeled by the geodesic 
motion of particles. From the geodesic motion on Finsler spacetimes we then derive 
the equations of motion for a Finsler fluid. 

The basic idea of the kinetic theory of fluids is the assumption that fluids are con¬ 
stituted by particles. In the simplest possible case of a single component fluid, which 
we will discuss here, all particles have identical properties, such as mass and electric 
charge, and follow piecewise geodesic curves. This geodesic motion corresponds to 
the motion of freely falling test masses without any other interaction. The inter¬ 
action between particles is modeled by collisions, which mark the endpoints of the 
geodesic pieces of the particle trajectories and correspond to instantaneous transfers 
of momentum between the particles. The physical background of this model of inter¬ 
actions is the assumption that the interaction distances are small compared to the 
distances between the particles and that the interaction times are short compared 
to the time between interactions, which is the case for a sufficiently low density. 

In order to construct a continuum theory of fluids from the particle model one 
introduces the one-particle distribution function </> : O —>■ such that for each 

oriented hypersurface a <Z O the integral 



(13) 


is the number of particle trajectories whose canonical lifts to O pass through a. 
Here particle trajectories 7 are counted positively (negatively) if the tangent vector 
of their canonical lift and a positively oriented basis of the tangent space to cr at 
the intersection point between 7 and a form a positively (negatively) oriented basis 
of the tangent space of O at that point. Further, uj denotes the particle measure 
introduced in the previous section. 

We now take a closer look at the canonical lifts of the particle trajectories to the 
observer space. It follows from our assumption of instantaneous momentum transfer 
that these lifts are discontinuous at collisions, i.e., they have endpoints correspond¬ 
ing to the particles’ velocities before and after the collision. To incorporate collisions 
into the continuum theory we define the collision transfer density </) : O —>■ K. such 
that for each hypervolume V C O the integral 



(14) 


is the number of initial points minus the number of final points of canonical lifts of 
particle trajectories in V. Here E denotes the volume form of the Sasaki metric G 
on O. 

Note that we have defined the counting prescription for curves passing a hyper¬ 
surface such that any curve which has an initial point, but no final point in the 
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hypervolume V passes its boundary dV in the positive direction, and in the nega- 
live direction in the opposite case. It thus follows that N[V] = N[dV]. From this 
we further find that 


<^S = N[V] = N[dV] = <l)Lo= , 


(15) 


IdV 


IV 


IV 


where we used Stokes’ theorem and the properties of the particle measure ui. Since 
this holds for any hypervolume V it follows that 


= Cr 


(16) 


The right hand side of this equation describes the evolution of (j). We thus obtain the 
equation of motion of the fluid by specifying a functional (j)[(j)] describing collisions 
between particles. For the simplest possible case of a collisionless fluid we have 
<p = 0, and the fluid equation of motion reduces to the Liouville equation Crip = 0. 

This concludes our discussion of the kinetic theory of fluids in general. In the 
following section we will discuss particular examples of fluids and derive their equa¬ 
tions of motion. 


5. Special cases 

After discussing general fluids on Finsler spacetimes in the previous section, we now 
consider a few special cases. The first case we display in section [5.II is dust, which 
is conventionally described as a collisionless perfect fluid with vanishing pressure. 
Starting from the Liouville equation we derive the generalization of the Euler equa¬ 
tions to dust on Finsler spacetimes. The second case shown in section 15.21 is the 
most general cosmological fluid. We start from a general Finsler spacetime with 
cosmological symmetry and derive the equations of motion for the most general 
fluid obeying the same symmetry. 


5.1. Dust fluid 


The first example we discuss is a collisionless dust fluid characterized by a matter 
density p{x) and four-velocity distribution u“(x), which is a future timelike vector 
field normalized by the Finsler function, and thus a function u : M ^ O. The aim of 
this section is to derive the equations of motion for p and u from the kinetic theory 
introduced in the preceding section. For this purpose we will first construct the one- 
particle distribution function p and then impose the Liouville equation Crp = 0 for 
a collisionless fluid. 

We start by introducing coordinates (x“,0“) = (x“,?/“/j/°) on O and writing 
the three-velocity as = u°‘ ju^. Together with the proper Dirac delta distribution 
on the unit timelike shell Sx we find the one-particle distribution function 


Q\ 1 - v{x)) 

p{x,9) = —p[x)—=== , 
m ^hP{x,e) 


( 17 ) 








August 14, 2015 0:22 WSPC/INSTRUCTION FILE archive'vl 


Non-metric fluid dynamics and cosmology on Finsler spacetimes 7 


where denotes the determinant of the restriction 

F F (.n. 

Q0a QQpSab 

of the Sasaki metric to S'^- A direct calculation of the Lie derivative £r</> is obstructed 
by the presence of the delta distribution. We therefore introduce an arbitrary smooth 
test function / : O —>■ R. Using the properties of the particle measure w and 
integrating by parts we then find 

0=/ d(0cc)/ = -/ dfAuj^=-[d{fL0)cj) = -[ 

Jo Jo Jo Jo Jo 

(19) 

Writing the Finsler function in the form 

F = y^Fix,y/y^)=y^Fix,e) (20) 

we obtain the Reeb vector field 


r = y 


da + F {N\e^ - iV“a) da 


( 21 ) 


using the partial derivatives da and da with respect to and and the volume 
form 


E = y/g^h^d^xd^e. 

Finally inserting the explicit form for </> we arrive at 

0 = [ d^xd^e^/Fsie - v)py<^ [4 + f {N\e^ - N^a) da 
Jo L 

Integrating over d^9 cancels the Dirac distribution and yields 

0= [ d‘^xWpufldaf + F{N\v‘^-N‘^a)daf] 


/• 


fM 


( 22 ) 

(23) 

(24) 




where it is indicated that all objects on observer space are taken at the point 
(x,v{x)). Observe that the partial derivative daf{x,v(x)) is fixed by f{x,v{x)) and 
thus not an independent quantity, while daf{x,v{x)) depends on the choice of / 
in a neighborhood of v{x), which can be chosen independently. We therefore need 
to eliminate daf using integration by parts. Since da is only a partial derivative, 
but f{x,v{x)) contains also an implicit dependence on x via v(x), we hrst need to 
rewrite the partial derivative into a total derivative using 


dafix,v{x)) = ^/(^)^ 4 )) - daV°'dafix,v{x)) . 


(25) 


Integration by parts and using the normalization u^F(x,v{x)) = F{x,u{x)) = 1 
then finally yields 


0 = / d‘^x\fF 

Jm 


N\—-N° 
1 


- 


Pdaf 


4(pw“) + (dagC + dau‘^ddg[a) 


(26) 


i/“=u“(4) 
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Since now all spacetime derivatives act on objects which depend only on spacetime 
coordinates, we can rename the coordinates back to and read off the equations 
of motion 

uUN\—-N‘^A-u^dau‘^ + ^-daU^ = 0 (27) 

and 

da{pu°-) + {dag[c + dau'^dd.gQ = 0 . (28) 

These equations can further be simplified. Using the properties of the Cartan non¬ 
linear connection we can rewrite the first equation (1271) in the form 

0 = u^{daU^ + N\) = Vu ^, (29) 

where we have introduced the dynamical covariant derivative V. Similarly, the sec¬ 
ond equation (1281) can be written as 

0 = 5a(pu“) -h ^pu°-g^’’‘'Sag[c = ^Sa{pu °‘), (30) 

where we have introduced the covariant derivative of the Cartan linear connec¬ 
tion. In the case of a metric Finsler function F'^{x, y) = \gab{x)y°‘y^\ these equations 
reduce to the well-known Euler equations 

= 0 and Va(pM“)=0 (31) 

for a pressureless fluid. 


5.2. Cosmological fluid 


As the hnal example we will derive the Liouville equation for a fluid with cosmologi¬ 
cal symmetry. For this purpose we introduce coordinates (t, r, d, (p) on the spacetime 
manifold M and the corresponding induced coordinates {t,r,'0,(p,y^,y'^,y'^,y'^) on 
TM. In these coordinates the cosmological symmetry generators take the form 


e /j- f ■ Q o , cosdcosip sin:^ \ 

\ r r sm v j 

^2 = ■\/l — kr^ ^sin d sin ipdr + 

^3 = \/l- kr"^ ^cosddr - —, (,6 = , 

^4 = sin ipd^ + > ^5 = - cos (pd^ + , (32) 

tan if tan if 

where k € { — 1,0,1} determines the spatial curvature of the corresponding space- 
time. For cosmologically symmetric fluid dynamics we require that both the back¬ 
ground geometry and the one-particle distribution function of the fluid obey this 
symmetry. 


cos d sin ip „ cos w _ 

- -dt/ + 

r r sm if 
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We start by deriving the most general cosmologically symmetric background 
geometry. A Finsler spacetime is symmetric under the action of a vector field ^ 
if and only if the geometry function L, and thus also the Finsle r function F, is 
invariant under the complete lift of ^ to the tangent bundle ,^2)3111] 

{Cda + y^dbCda) F = 0. (33) 

In the present case of a cosmological symmetry it is most convenient to introduce 
new coordinates on TM via the definition 

t = i, r = r, ■d = fl, = = y, 

y'" = w cos uyj 1 — = — sin u cos u , y^ = -^ sin u sin D . (34) 

f f sin i} 

After calculating the complete lifts of the vector fields (l32ll in these coordinates 
it turns out that the most general Finsler function with cosmological symmetry is 
given hy F = F{t,y,w). Further, from the fact that F is 1-homogeneous in the 
coordinates y and w follows that it must be of the form 

F(f, y, w) = yP (£, w/y) . (35) 

For the discussion of geodesic motion it is convenient to introduce yet another set 
of coordinates on TM as 

t = t, f = f , i} = i} , ifi = (f, u = it, V = V , y = yF t, — , w = — . 

\ y J y 

(36) 

In these coordinates the observer space O is given as a connected component of 
the submanifold y = 1, so that one can use the remaining seven coordinates to 
parametrize O. The Reeb vector field then takes the form 


= — ( + ih cos u\/l — kF^dr + 


w sm u cos V 


F 




w sm u sm v . 


w sin 


du - 


fsinfl 
w sin u sin v 


'land 


p: ^tw a 

Uy Oy 


, (37) 


where the subscripts t and w indicate derivatives with respect to t and w, respec¬ 
tively. 

We finally come to the discussion of fluids on the Finsler spacetime background 
derived above. Since the background geometry obeys the cosmological symmetry 
defined by the vector fields ([32l) . their canonical lifts are tangent to the observer 
space O. A fluid obeys the same symmetry if and only if its one-particle distribu¬ 
tion function (j) is invariant under the restriction of these canonical lifts to O. In 
the present case the most general one-particle distribution function satisfying this 
condition takes the form ip = (Pitjw). Its Lie derivative with respect to the Reeb 
vector field is thus given by 

Cr(j)= ^ ((pt-■ (38) 

^ V ^WW J 
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For the simplest possible case of a collisionless fluid the equation of motion hence 
takes the form 

^ww4^t - ^tw4*w ■ (^ 9 ) 

This is the Liouville equation for a fluid with cosmological symmetry. 

6. Conclusion 

We have derived a model for fluid dynamics on Finsler geometric backgrounds based 
on the kinetic theory of fluids and applied our model to two important special cases. 
Our results show that Finsler spacetimes provide a suitable background geometry 
for fluid dynamics and that the obtained fluid equations of motion reduce to their 
well-known limits if the background geometry is metric. 

Of course any model of fluid dynamics must be complemented by a suitable 
model of gravitational dynamics in order to derive consistent solutions. I^or the 
Finsler spacetimes we used here an extension of general relativity exists.^^El xhe 
source of the gravitational field in this model is a scalar function on the tangent 
bundle of spacetime. Deriving this energy-momentum scalar for a kinetic fluid and 
the corresponding gravitational field equations will be a topic of future research. 
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